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INTENSE  ION  BEAM  GENERATION, 
PLASMA  RADIATION  SOURCE  AND 
PLASMA  OPENING  SWITCH  RESEARCH 


D.A.  Hammer,  M.D.  Coleman,  N.  Qi,  P.L.  Cimilon,  and  R.N.  Sudan 


Final  Report  -  ONR  Contract  N00014-88-K-2009 


ABSTRACT 

This  report  describes  research  on  intense  ion  beam  diodes,  plasma  opening 

> 

switches  and  dense  z-pinch  plasma  radiators.  Laser-induced  fluorescence  spec¬ 
troscopy  has  been  used  to  map  the  electrostatic  potential  profile  in  a  plasma^prefilled 
magnetically  insulated  ion  diode.  In  a  simple  planar  diode,  the  measured  profile  is 
inconsistent  with  the  electrons  being  confined  in  a  sheath  near  the  cathode  by  the 
magnetic  field.  Rather,  the  profile  implies  the  presence  of  electrons  throughout  the 
accelerating  gap.  A  theoretical  model  of  the  penetration  of  current  and  magnetic 
field  into  a  plasma,  and  of  the  current-driven  effective  collision  frequency  has  been 
developed.  The  snowplow  action  of  the  rising  magnetic  field  causes  a  steep  rise 
in  the  plasma  density  at  the  leading  edge.  The  subsequent  multistreaming  of  the 
ions  caused  by  ion  reflection  at  the  current  layer  could  lead  to  ion  heating  through 

l 

collective  effects.  The  two-dimensional  electron  flow  in  the  plasma-cathode  vacuum 
gap  is  also  treated.  Dense  z-pinch  plasma  radiation  source  experiments  have  been 
initiated  on  the  LION  accelerator  using  gas  puff  and  fine  wire  loads.  The  x-pinch 
was  found  to  be  a  more  effective  way  to  generate  soft  x-rays  than  a  single  wire 
pinch  or  a  gas  puff  implosion.  Plasma  opening  switch  experiments  being  initiated,  — 
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and  plasma  anode  ion  diode  development  work  being  terminated  are  also  briefly 
described. 

I.  INTRODUCTION 
I.A.  Summary 

During  the  period  1  January  1988  to  31  December  1988,  research  was  in 
progress  at  Cornell  University  under  the  support  of  the  Naval  Research  Labora¬ 
tory  (ONR  Contract  N00014-88-K-2009)  on  five  research  tasks.  They  were: 

1.  Design  and  test  an  active  anode  plasma  ion  source  for  a  500  kV  magnetically 
insulated  diode  to  be  operated  on  the  3fl,  120  ns  Neptune  Pulsed  Power  Gen¬ 
erator; 

2.  Carry  out  spectroscopic  investigations  of  conditions  in  detail  in  the  high  voltage 
gap  of  a  magnetically  insulated  ion  diode; 

3.  Initiate  experiments  on  a  long  conduction  time  plasma  opening  switch; 

4.  Initiate  theoretical  investigations  of  plasma  opening  switches; 

5.  Initiate  experiments  and  calculations  on  z-pinch  x-ray  sources  with  particular 
application  to  photo-pumped  or  recombination-pumped  XUV  lasers. 

In  this  report  we  briefly  describe  our  progress  on  items  1  and  3  since  the  first  is 
still  in  progress,  but  is  no  longer  supported  by  NRL,  while  the  third  is  just  beginning. 
These  two  are  discussed  in  this  introduction.  More  complete  presentations  of  results 
on  tasks  2,  4,  and  5  are  presented  in  Sections  II,  III,  and  IV,  respectively.  Tasks  4 
and  5  were  also  partly  supported  by  the  DOE  Office  of  Inertial  Fusion. 

I.B.  Active  Anode  Plasma  Ion  Source  for  Neptune 

The  motivation  for  this  ion  source  developmend  task  was  our  successful  oper¬ 
ation  of  an  active  anode  plasma  ion  source  on  the  150  kV,  1/zs  pulse,  lfl  LONG- 
SHOT  pulsed  power  generator  under  NRL  support.  The  completion  of  this  task 
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will  demonstrate  both  higher  voltage  and  shorter  pulse  operation  of  such  a  “plasma 
anode  diode.”  As  of  the  end  of  the  period  covered  by  this  final  report,  the  design 
and  construction  of  the  hardware  for  the  plasma  anode  diode  was  completed,  and 
vacuum  testing  was  carried  out.  The  diode  contains  five  independent  high  voltage 
electrical  circuits  which  must  operate  correctly  to  produce  the  plasma  in  the  de¬ 
sired  diode  magnetic  insulation  field  configuration,  and  electrical  testing  of  these 
circuits  had  begun.  The  work  is  continuing  under  DOE  sponsorship  since,  by  mu¬ 
tual  agreement  of  the  Cornell  investigators  and  NRL  technical  personnel,  this  work 
is  no  longer  sponsored  by  NRL. 

I.C.  Long  Conduction  Time  Plasma  Opening  Switch 

The  success  of  relatively  small  and  inexpensive  pulsed  power  systems  based 
on  inductive  energy  storage  depends  on  the  development  of  efficient  fast  opening 
switches.  The  long  conduction  time  plasma  opening  switch  (POS)  is  a  good  candi¬ 
date  for  the  last  stage  opening  switch  in  such  a  system.  A  basic  understanding  of 
the  mechanisms  involved  in  opening  are  needed  to  optimize  the  behavior  of  these 
switches.  At  present  they  are  characterized  by  longer  switching  times  and  more 
losses  than  the  short  conduction  time  switches. 

We  have  recently  begun  experiments  similar  to  earlier  work  at  NRL  (D.D. 
Hinselwood  et  al.,  Appl.  Phys.  Lett.  49,  1635  (1986))  to  investigate  the  physics 
of  the  operation  of  long  conduction  time  opening  switches.  Our  emphasis  will  be 
on  a  determination  of  conditions  in  detail  using  spectroscopic  diagnostics,  although 
we  will  also  use  the  usual  magnetic  probes  to  monitor  current  penetration  and 
Faraday  cups  for  plasma  flux.  The  complete  system  was  assembled  prior  to  the  end 
of  1988  and  initial  experiments  were  carried  out.  The  results  are  summarized  in 
the  following  paragraph.  This  research  is  continuing  under  a  new  NRL-supported 
grant. 
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The  POS  system  is  pulsed  by  a  1.9/xF  Scyllac  capacitor  charged  to  50  kV.  It 
has  an  inductance  of  300  nH  (no  plasma  in  the  switch,  short  circuit  load)  with  10 
nH  of  this  downstream  of  the  POS.  The  plasma  for  the  10  cm  x  10  cm  planar 
switch  is  provided  by  a  Mendel-type  carbon  plasma  gun  which  is  pulsed  by  a  2.0#iF 
capacitor  charged  to  20  kV.  This  supplies  a  plasma  with  roughly  3  xlO14  cm'3 
ion  density,  predominantly  C++ ,  with  some  C+++,  with  an  electron  temperature 
of  10-15  eV.  In  a  preliminary  run  with  this  switch  plasma,  conduction  times  up 
to  400  ns  were  observed  in  a  short  circuit  load.  An  opening  time  of  130  ns  was 
obtained,  delivering  70  kA  (out  of  106  kA  available).  An  electron  beam  diode  load 
gave  220  ns  conduction  time  and  190  ns  opening  time  with  70  kA  (out  of  100  kA 
available)  delivered  to  the  diode  with  some  voltage  multiplication.  The  dynamics 
of  POS  operation  will  be  studied  in  the  next  experimental  series. 

II.  DETERMINATION  OF  CONDITIONS  IN  DETAIL  IN  A  PLASMA- 

PREFILLED  ION  DIODE  BY  SPECTROSCOPIC  TECHNIQUES 

H.A  Introduction 

Emission  spectroscopy  and  Laser  Induced  Fluorescence  spectroscopy  (LIF)  pro¬ 
vide  potentially  valuable  techniques  for  measuring  conditions  in  high  voltage  gaps 
non-intrusively.  In  previous  work,  the  electric  field  was  measured  in  a  magnetically 
insulated  ion  diode  (MID)  using  Stark  shift  emission  spectroscopy1.  In  this  paper, 
we  report  measurements  of  the  electrostatic  potential  profile  in  a  plasma-prefilled 
MID  using  LIF.  An  accelerating  gap  of  3.5  mm  is  shown  to  form  in  1012-1014  cm-3 
argon  plasma  in  a  few  ns  when  a  350  kV  pulse  is  applied  by  a  lOfJ  pulsed  power 
generator.  The  observed  potential  profile  for  a  well  insulated  gap  (B/Bcr,t  =  1.75, 
where  B  is  the  applied  magnetic  field  and  B crit  is  the  minimum  field  for  magnetic 
insulation2  of  the  actual  accelerating  gap)  is  inconsistent  with  the  electrons  being 
confined  in  a  thin  sheath  by  the  insulating  magnetic  field. 
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n.B  The  Plasma-Prefilled  Ion  Diode 

Figure  ILl  shows  an  end  view  of  the  MID.  The  anode  is  connected  to  the  center 
conductor  of  a  10  0  pulse  line.  Argon  gas  is  puffed  into  the  interior  of  the  anode. 
It  is  ionized  by  the  electric  field  induced  by  the  rapidly  rising  insulating  magnetic 
field  of  the  MID.  The  magnetic  field  is  generated  by  current  pulsed  through  the 
cathode,  which  acts  as  a  single  turn  coil.  The  plasma  passes  through  slots  in  the 
anode  (oriented  vertically  in  Fig.  II. l)  and  then  into  the  accelerating  gap.  The 
slotted  area  is  6  cm  high  by  10  cm  long.  A  variable  delay  adjusts  the  time  of  arrival 
of  a  400  kV,  80  ns  pulse  from  the  10  Cl  pulsed  power  generator  relative  to  when  the 
plasma  is  formed.  This  allows  the  peak  plasma  density  and  the  density  profile  to  be 
varied  as  part  of  this  experiment.  Plasma  density  in  the  acceleration  gap  typically 
ranges  from  —  1014  cm-3  near  the  anode  to  ~  1012  cm-3  near  the  cathode  at 
the  time  the  high  voltage  pulse  arrives.  Density  is  determined  from  the  absolute 
intensity  of  the  4806A  Ar+  line.  The  spectrometer /photomultiplier  system  used  for 
this  measurement  was  calibrated  by  simultaneously  measuring  the  plasma  density 
with  a  Langmuir  probe  at  one  particular  location.  The  plasma  temperature  was 
3.5  eV,  as  measured  from  Ar+  line  intensity  ratios. 

An  argon  plasma  is  used  for  these  experiments  because  of  the  existence  of  a 
metastable  state  (3d4D)  in  Ar+  16.4  eV  above  the  ground  state.  Pumping  from 
a  ground  state  or  a  metastable  state  is  necessary  for  LIF  studies  at  these  plasma 
temperatures  and  densities,  and  pumping  from  the  ground  state  is  difficult  due  to 
the  large  energy  separation  to  the  first  excited  state  in  most  ions. 

Figure  II.2  shows  diode  voltage,  current,  and  calculated  impedance  time  histo¬ 
ries  for  a  shot  with  a  relatively  high  density  initial  plasma-prefill  (about  1013  cm-3 
near  the  cathode).  There  is  an  initial  period  of  low  diode  impedance  before  the 
voltage  begins  to  rise.  By  firing  the  diode  with  less  plasma  in  the  acceleration  gap, 
significant  variation  in  the  impedance  characteristics  can  be  obtained,  including 


5 


Figure  II.l  End  view  of  the  plasma-prefilled  ion  diode. 

cases  in  which  the  diode  impedance  is  relatively  constant  throughout  the  voltage 
pulse.  Typically  30-60  A/cm2  Ar+  beams  axe  extracted  as  determined  from  biased 
Faraday  cups. 

II. C  LIF  Electrostatic  Potential  Measurement 


The  experimental  setup  is  illustrated  in  Fig.  II.3.  A  dye  laser,  tuned  near  the 
440lA  Ar+  line,  is  synchronized  with  the  high  voltage  pulse  and  directed  through 
slots  in  the  MID  cathode  (also  vertical  in  Fig.  II.l)  toward  the  anode.  The  laser 
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Figure  U.2  Time  dependence  of  (from  top  to  bottom)  diode  current,  voltage,  and 
calculated  impedance,  with  scales  for  each  indicated  on  the  ordinate. 

light  excites  electrons  from  the  3d4D  metastable  level  of  the  ions  to  the  4p4P  level, 
from  which  the  ions  fluoresce  at  4806A.  The  fluorescence  emitted  perpendicular  to 
the  ion  beam  acceleration  direction  and  along  the  magnetic  field  lines  is  directed 
into  a  spectrometer  and  monitored  by  a  photomultiplier  tube.  The  wavelength 
for  absorption  of  the  laser  light  by  the  metastable  Ar+  ions  is  a  function  of  the 
distance  from  the  anode,  x,  since  the  increasing  ion  velocity,  v,  Doppler  shifts  the 
wavelength  for  resonant  absorption  to  values  which  increase  with  x.  To  determine 
the  ion  velocity  at  a  particular  x,  a  series  of  shots  is  taken  tuning  the  laser  to  longer 
and  longer  wavelength  above  440lA.  The  detuning  from  440lA,  AA,  above  which 
no  fluorescence  is  detected,  then  determines  the  velocity  at  that  x: 

v  =  (AA/440l)c. 

The  potential  drop  relative  to  the  anode  potential  is  th<  n  given  by  $(x)  =  |mu 2, 
where  m  is  the  mass  of  an  argon  atom,  and  c  is  the  speed  of  light.  A  complete 
$(x)  profile  is  obtained  by  repeating  the  procedure  with  the  collection  optics  for 


Figure  n.3  Schematic  diagram  of  the  experimental  arrangement 

the  spectrometer  focussed  at  different  values  of  x. 

In  the  following  discussion  of  results,  the  laser  beam  was  focussed  using  a  cylin¬ 
drical  lens  so  that  it  was  about  6  cm  wide  and  about  3  mm  high  in  the  acceleration 
gap.  The  15  ns  wide  (FWHM)  pulses  contained  12  mJ.  The  x  resolution  of  the 
collection  optics  was  0.5  mm.  Figure  II.4  shows  a  typical  plot  of  fluorescence  vs. 
detuning  of  the  laser  above  the  unshifted  wavelength  for  a  single  x  position.  Notice 
that  the  laser-induced  fluorescence  drops  to  zero  in  1-1. 5  A  around  A  A  =  13. 5A. 
This  gives  v  =  9.2  x  107  cm/s  and  $  =  176  kV  for  this  particular  x. 

Figure  II.5  shows  a  time  history  of  the  potential  profile  for  shots  with  voltage 
and  current  waveforms  like  those  in  Fig.  II. 2,  and  an  insulating  B-field  of  7  kG. 
The  times  10,  23,  37,  and  50  ns  are  referenced  to  the  time  that  the  current  pulse 
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12 


Laser  detuning  (A) 

Figure  II.4  Typical  plot  of  fluorescence  vs.  laser  detuning. 

starts  to  rise  in  the  diode.  (The  timing  uncertainty  is  about  2  ns.)  Notice  that 
there  is  a  >  8  mm  wide  region  adjacent  to  the  anode  which  is  at  anode  potential. 
All  of  the  ion  acceleration  occurs  in  a  3.5  mm  wide  region  adjacent  to  the  cathode. 
The  points  at  the  cathode  are  the  inductively  corrected  voltage  obtained  from  a 
capacitive  probe  on  the  end  of  the  pulseline.  Such  voltage  measurements  are  always 
consistent  with  the  ion  kinetic  energies  measured  by  the  LIF  technique  after  the 
ions  have  been  extracted  from  the  diode  within  10%. 

II.D  Discussions  and  Conclusions 
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Figure  13.5  Time  dependence  of  ion  kinetic  energy  (which  equals  the  potential  drop 
from  the  anode)  as  a  function  of  distance  from  the  anode.  At  t  =  50  ns,  B/Bcr,t  = 
1.1.  The  metal  anode  is  at  x  =  0. 

The  data  shown  in  Fig.  H.S  was  selected  for  presentation  because  the  density 
was  high  enough  that  the  opening  up  of  the  acceleration  gap  was  visible.  The 
corresponding  voltage  waveform  is  shown  in  Fig.  II.2.  Thus,  the  impedance  history 
of  the  gap  during  the  first  15  ns  is  similar  to  that  observed  for  a  plasma  opening 
switch. 

Potential  profiles  have  also  been  obtained  for  a  plasma-prefill  about  a  factor  of 
three  lower  density  than  was  the  case  when  the  profiles  in  Fig.  II.5  were  obtained. 
In  that  case,  the  diode  voltage  and  current  rise  together.  Finally,  we  have  also 
investigated  diode  configurations  in  which  electron  emitting  vanes  have  been  added 
to  the  cathode.  In  cases  in  which  the  effective  diode  gap  (the  gap  between  the 
cathode  and  *he  largest  x  at  anode  potential)  is  well  insulated  (e.g.  -  B/Bcrit  = 
1.75),  the  potential  profile  clearly  does  not  correspond  well  to  models  involving 
a  well  confined  electron  sheath.3  Instead,  the  profile  corresponds  nicely  with  the 
simple  assumption  of  uniform  electron  density  in  the  accelerating  gap,  in  qualitative 
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agreement  with  the  electric  field  measurements  made  in  a  surface  flashover  MID 
using  emission  spectroscopy  by  Maron  et  al.1 

This  LIF  potential  measuring  technique  could  possibly  be  utilized  in  other, 
higher  power  ion  diodes  which  cannot  use  argon  gas  as  the  anode  plasma  source,  by 
seeding  the  plasma  with  desirable  impurites  (e.g.-  calcium,  in  which  pumping  from 
the  ground  state  of  the  ion  is  possible).  Also,  this  technique  could  possibly  be  put 
to  great  use  as  a  non-intrusive  method  to  measure  voltages  (as  small  as  2  kV)  as  a 
function  of  all  spatial  directions  in  plasma  opening  switches. 
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m.  CURRENT  PENETRATION  AND  ION  DYNAMICS  IN  A 

PLASMA  OPENING  SWITCH 

m.A  Introduction 

Progress  in  the  physics  and  technology  of  plasma  opening  switches  (POS)  for 
pulsed  power  systems  in  the  terawatt,  less  than  100  nanosecond,  range  hats  been 
reported  recently  in  the  Dec.  1987  issue  of  Plasma  Science.1  There  is,  however, 
considerable  controversy  over  the  physical  mechanisms  underlying  the  operation  of 
this  switch.  In  the  main  there  axe  two  interpretations  of  switch  operation.  The  NRL 
group2  maintains  that  in  the  later  stages  the  electrons  are  magnetically  insulated 
by  their  self-magnetic  field  and  are  unable  to  cross  the  gap  between  the  cathode  and 
the  plasma.  The  conduction  current  is  therefore  provided  by  the  plasma  ions  and 
because  this  is  greatly  in  excess  of  the  ion  thermal  flux,  results  in  the  erosion  of  the 
plasma  surface  facing  the  cathode  thereby  increasing  the  A-K  gap  and  eventually 
to  current  interruption.  The  alternative  picture3  argues  that  the  magnetic  j  x  B 
forces  must  play  a  strong  role  in  the  current  interruption  process. 

In  this  paper  we  address  the  issues  of  (i)  current  penetration  into  a  plasma 
from  a  high  power  travelling  electromagnetic  pulse  incident  on  the  surface  of  plasma 
created  in  a  transmission  line,  (ii)  the  ion  dynamics,  and  (iii)  the  two-dimensional 
electron  flew  in  the  plasma/cathode  vacuum  gap.  A  proper  understanding  of  these 
problems  is  an  absolute  requisite  for  any  POS  model. 

III.B  Mathematical  Model 

We  consider  a  one-dimensional  representation  of  the  plasma,  deferring  the  con¬ 
sideration  of  cathode  and  anode  sheaths  to  Section  VI.  Plasma  quasineutrality  is 
assumed  i.e.,  ne  =  =  n;  the  ions  are  taken  to  be  singly  ionized  without  much 

loss  in  generality.  In  the  two- fluid  approximation  the  equations  governing  plasma 
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behaviour  are  given  by: 


Here  p  =  nM  is  the  plasma  mass  density,  vt  =  =  v  is  the  x-component  of  the 

ion  and  electron  velocities,  w  is  the  electron  y- velocity,  p  —  mqtv  —  eA/c  is  the 
electron  y-canonical  momentum,  7  =  (l  —  (v2  +  w2)/c2)  ^  ,  B  =  B  2  and  the 
vector  potential  A  =  A  y.  The  momentum  transfer  between  electrons  and  ions  is 
determined  by  the  collision  frequency  u  which  signifies  both  Coulomb  collisions  and 
collective  effects  due  to  microinstabilities.  Equation  (3)  expresses  the  conservation 
of  the  electron  y-canonical  momentum  if  the  collisions  are  infrequent.  The  pressure 
is  neglected  because  the  initial  plasma  /?  in  the  current  sheath  is  small. 

In  the  collisionless  limit  dp/dt  »  mu  tv,  Eqns.  (3)  and  (4)  may  be  combined 
to  give 


d2A 
dx 2 


4n  ne2A 
me 2 


(5) 


In  the  non-relativistic  limit  eA  «  me2  this  equation  determines  the  current  sheath 
thickness  A  to  be  of  order  \e  —  c/we,  u2  =  4x n0t2  jm.  In  the  ultra-relativistic 
limit  eA  »  me 2  and  B  =  Bq(x  =  0)  —  47rn0ea:.  This  solution  is  limited  to  t  < 
A/cx,  [ca  =  Bo/(4irpo)3  is  the  Alfv4n  velocity],  during  which  the  density  n  does 
not  change  significantly  from  the  initial  value  no.  Note  that  in  the  ultrarelativistic 
limit  A  =  Bo/4nn0e. 

In  almost  all  instances  of  pulsed  power  POS  applications  the  electron  drift  ve¬ 
locity  in  the  current  sheath  far  exceeds  the  acoustic  speed  c,  —  ( T / M )  $  and  perhaps 
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even  approaches  or  exceeds  (T/m)^  the  electron  thermal  speed.  The  microinstabil- 
ities  excited  by  the  current  lead  to  an  enhanced  rate  of  momentum  transfer  from  the 
electrons  to  the  ions  or  to  the  emitted  waves.4  This  increase  in  the  effective  collision 
frequency  t/*  may  cause  the  collisional  term  in  Eqn.  (3)  to  dominate  the  electron 
mechanical  inertia.  If  this  limit  obtains  Eqns.  (3)  and  (4)  may  be  combined  to  give 


dB  d_ 
dt  dx 


3J3 
dx ’ 


(6) 


where  rf  —  (c2rm/*/47rne2).  Thus  the  collisionless  limit  is  defined  by  Eqns.  (1),  (2) 
and  (5)  and  the  collisional  limit  by  Eqns.  (1),  (2)  and  (6). 

ULC  Microinstability  in  the  Current  Sheath 

A  recent  calculation  by  Kulsrud  et  al.s  explores  the  consequences  of  ion  acoustic 
waves  driven  unstable  by  the  electron  drift.  They  show  that  wave  emission  by  the 
electrons  leads  to  a  drag  force  on  the  electrons  in  the  y-direction  which  is  equivalent 
to  an  effective  collision  frequency 


u.  =  (*/2)*  (*;/*„)  (T/m)  *  («  v/T)\  (7) 


where  the  spectrum  has  a  mean  wave- vector  k  =  (0,fcy,fc||)  and  p  is  the  amplitude 
of  the  electric  potential  fluctuations. 

If,  as  suggested  by  Kulsrud  et  a/.5,  we  set  ky ~  1,  (A©  is  the  Debye  length) 
fc||t»e  ~  ky w,  ( ve  and  w  are  the  electron  thermal  and  drift  velocities),  and  ejt/T  ~  1, 
one  obtains  i/«  ~  we(ve/u»),  where  we  is  the  electron  plasma  frequency.  At  the 
high  saturation  amplitude  chosen  t4>jT  ~  1  the  nonlinear  distortion  of  the  electron 
orbits6  by  the  emitted  wave  fields,  omitted  in  the  earlier  analysis,5  must  be  taken 
into  account.  In  the  Appendix  we  employ  the  nonlinear  drift-kinetic  equation  for 
the  electrons  and  compute  the  renormalization  of  the  wave  growth  rate  due  to 
stochastic  diffusion  of  electron  trajectories  in  the  emitted  fields.  We  also  consider 
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the  nonlinear  exchange  of  momentum  between  the  waves  and  the  ions,  which  is 
expressed  in  terms  of  an  anomalous  “eddy”  viscosity.  At  the  saturation  intensity, 
the  ioii.  damping  balances  the  electron  excitation,  and  therefore  limits  the  fluctuation 
amplitude  to 

c4>/T  ~  {Clelkve)x/*(w/ve)l/3  <  0(1);  1  >  kve/Ue  »  w/ve.  (8a) 


The  momentum  transfer  rate  from  electrons  to  waves  is  reduced  due  to  nonlinear 
scattering  with  respect  to  the  linear  estimate.  Altogether,  it  is  found  that  the 
effective  collision  frequency  is  a  factor 

i(n,/ivt)4'3(w/v,)‘/3  «  1  (86) 

smaller  than  Kulsrud  et  al.’s5  estimate. 


D3.D  Solution  of  Model  Equations 


Equations  (1),  (2),  (5),  and  (6)  may  be  expressed  in  Lagrangian  variables 
through  a  coordinate  transformation  (x,  t)  <=>  (£(x,t),  r)  such  that  £  =  x  at  t  = 
0  and, 


dt 

dt 


H 

dt 


+  vTxi  =  a' 


(9) 


r  =  t. 


(10) 


Thus,  because  of  density  conservation, 


dijdx  =  n(£,r)/n0(£)> 


(11) 


where  n0(£)  is  the  initial  density  profile  at  r  =  0.  In  the  new  coordinates  Eqns.  (l) 
and  (2)  become 


d  .  .  .  dv 

Tt  (no/n)  "  T(  =  °- 

dv  1 _  d_  j&_  _ 

dr  no  d£  8 irM 


(12) 

(13) 
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Combining  (12)  and  (13)  we  get 


d 2  /1\  l_d_  l_  d_  B2 

dr2  \n/  no  no  8irM 

For  the  collisionless  case  Eqn.  (5)  transforms  to 

dB  _  47rn0(Oc2 
d£  me 2 

In  the  ultrarelativistic  case  the  exact  solution  of  (15)  is 

J9(£)  =  Bq  —  4ire  f  n0d£ 


(14) 


(15) 


(16) 


even  when  the  plasma  dynamics  is  taken  into  account.  For  the  collisional  case,  Eqn. 
(6)  transforms  to 

57  =  §1  *°  5?’ ■  (17) 

where  rio  —  c2mv*/4irn0e2  and  B  at  £  =  0  is  equal  to  Bq(t),  a  given  function 
determined  by  the  imposed  current  pulse.  In  terms  of  the  vector  potential  A,  Eqn. 
(17)  may  be  rewritten  as 


dA 

dr 


(18) 


Ultra-relativistic  Collisionless  Case 

We  treat  this  case  because  of  its  simplicity,  its  relevance  for  very  high  currents 
and  also  because  it  reveals  a  very  important  feature  of  this  problem.  Substituting 
Eqn.  (16)  in  Eqn.  (14)  yields, 


<P_ 

dr2 


4?re2 

+  -aT  “  °- 


(19a) 


thus 


n(£,  t)  =  »0(£) 


1 


-  w?  r2 
2  * 


-l 


(196) 


16 


for  £  <  £c,  where  Bo  =  4 ire  f£a  n0d(  and  ur?  =  4irno(£)c2  /M.  We  notice 
immediately  that  a  singularity  arises  in  the  density  at  rj  =  V^/W  This  singularity 
signals  the  breakdown  of  the  fluid  description  and  this  occurs  in  a  time  much  less 
than  a  plasma  ion  period!  The  flow  velocity  is  given  by 

t/(£,  r)  =  eBo/M  -  ( 4tcc2/M )  J  d£  no  r.  (20) 

Collisional  Case 

The  singularity  in  density  arises  also  in  the  collisional  case.  From  Eqn.  (17) 
with  B(£  =  0,  r)  determined  by  the  external  current  as  a  monotonically  increasing 
function  of  r,  it  is  apparent  that  J9(£,  r)  will  be  a  monotonically  decreasing  function 
of  £  for  fixed  r,  because  of  the  diffusive  nature  of  this  equation.  Consequently  the 
second  term  of  Eqn.  (14)  will  be  positive  definite.  Thus 

d2  ..  ,  1  d  1  d  B2  , 

dr 2  n0  n0  d£  8 ~  K^'  ^ 

and 

T  T’  -1 

»(£,  r)  =  n0(£)  1-/  dr'f  dr"  K{t,r")  .  (21) 

Jo  Jo 

From  dimensional  consideration  of  Eqns.  (14)  and  (17)  the  approximate  order  of 
time  needed  for  the  singularity  to  develop  is  r2  «  2«/./nefi,,  where  fle  and  ft,  are 
the  electron  and  ion  cyclotron  frequencies  in  Bo-  In  the  non-relativistic  collisionless 
case  rj  «  2jy/Cie^li.  Thus,  if  r2  <  n,  i.e., 

v./2y/n eft,  <  1,  (22) 

the  situation  is  treated  as  collisionless  and  collisional  if  the  above  inequality  is 
violated. 

The  physical  implication  of  the  singularity  in  density  is  that  the  snowplowing 
action  of  the  magnetic  field  results  in  the  generation  of  a  reflected  stream  of  ions 
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from  the  current  layer  because  of  the  electrostatic  potential  in  the  current  layer 
given  by 


dip  _  w  dA  _  1  d 
d£  c  d£  87rn0e  d£ 

An  approximate  expression  for  the  potential  drop  in  x  is  epo 
ultrarelativistic  limit  p  «  A). 


(23) 


Bq /8tt n0e  (in  the 


m.E  Vlasov  Approach 


The  effect  of  thermal  spread  in  the  ion  dynamics  is  obtained  through  the  Vlasov 

equation  for  the  ion  distribution  f(vx,  vyt  x,  t)  which  for  our  problem  becomes 

df  df  e  /  dp  vy  dA\  df  t  / dA  eE4\  df  _ 

dt  Vx  dx  +  M  \  dx  c  dx )  dvx  Me  \  dt  +  Vx  dx )  dvy 

(24) 

We  may  neglect  (vy/c)  dA/dx  as  compared  to  —  dpjdx ,  (see  Eqn.  (23))  and  employ 
the  transformation  to  the  Lagrangian  coordinate  £  as  in  Eqn.  (9).  Expressing  / 
also  as  a  function  of  the  ion  canonical  y-momentum  P  =  Mvy  +  eA/c,  Eqn.  (24) 
reduces  to 


df  .df  Id  tB\  df 

dr  Vx  V  d£  n0M  d£  8?r  dvx 

The  characteristics  of  this  equation  are  determined  by 


0. 


(25) 


nod£/n[vx  —  v(^,  r)]  =  dvx/(d/nod£)B2 /8nM  =  dr, 


from  which  we  obtain 

t  +  +  2«  Tr  'n(”°/n)  =  "^57  +  ^  (26) 

In  the  hydrodynamic  limit  with  no  thermal  spread,  the  ions  do  not  move  with 


respect  to  the  background  plasma,  £  —  0  =  £,  and  we  recover  Eqn.  (13).  The 


trajectory  of  a  test  particle  which  has  a  thermal  velocity  Svx  in  excess  of  the  hy¬ 


drodynamic  velocity  v{£,r)  at  r  =  0  is  obtained  from  Eqn.  (26)  to  be 


£  =  6vx{n/n0 )  exp  - 


r 


(27) 
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where  n/fio  is  a  function  of  $,  r  obtained  from  the  hydrodynamic  equations.  In  the 
ultraxelativistic  limit,  (n/n0)  is  not  explicitly  dependent  on  £,  Eqn.  (19b)  provided 
0  <  £  <  £c,  and  Eqn.  (27)  can  be  readily  integrated  to  give  £  =  Svx(-^)2.  Thus 
the  ion  temperature  ~<  (£«o/n)2  >  increases  dramatically. 

After  r  =  T\  or  r2  the  time  taken  for  a  singularity  to  develop  in  the  density, 
multistreaming  takes  place  and  our  Lagrangian  transformation  which  depends  on 
a  one  to  one  correspondence  between  x  and  £  is  no  longer  valid.  The  effect  of  a 
thermal  spread  in  the  ions  is  to  reduce  the  magnitude  of  the  density  increase  because 
all  the  particles  will  not  arrive  at  the  same  time  at  the  same  point. 

The  physical  consequences  of  ion  reflection  at  the  current  layer  after  r  = 
Ti  or  r2,  is  the  development  of  ion  viscosity  between  the  incoming  and  reflected 
ions  as  a  result  of  micro-instabilities  such  as  the  lower-hybrid.  However,  these  in¬ 
stabilities  will  only  occur  in  the  current  layer  where  a  finite  magnetic  field  exists.  If 
the  ion  slowing  down  length  (as  a  result  of  ion  viscosity)  is  greater  than  the  current 
layer  width  A,  the  reflected  ions  will  proceed  in  ballistic  orbits  in  the  magnetic 
field-free  region  of  the  plasma. 

Because  of  the  complications  introduced  by  multistreaming  after  r  =  riorr2, 
analytical  solutions  are  abandoned  in  favor  of  a  numerical  solution.  The  basic 
equations  for  a  simple  one-dimensional  hybrid  code  under  development  in  which 
the  electrons  are  treated  as  a  massless  fluid  and  the  ions  as  discrete  particles,  are 
as  follows: 


ij'  ~  Me  {(Ui 


+ 


dA 


4ir  ne  dx 2  /  dx 


+ 


.  3A) 


(28) 


e  :  e  (dA 

Vj  ~  ~  WcA  +  Wc\~di  + 


e  / dA  dA\ 

4c\dt  U*  dx ) 


~  Me  ~  Ux)  dP  (29) 


dA  mc^v  d2A 
dx  4tt  ne2  dx 2 


dA 

dt 


+ 
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(30) 


where  the  mean  velocities  of  particle  ions 


C:)  -  i  t  g) 


at  the  location  x  and  n  is  the  mean  ion  density. 


m.F  Two-Dimensional  Electron  Flow  in  the  Plasma/Cathode 
Vacuum  Gap 

We  show  below  that  the  current  between  the  cathode  and  the  plasma  in  the 
vacuum  gap  can  be  carried  by  the  electrons  even  when  strongly  magnetized  in  a 
two-dimensional  equilibrium  flow.  We  have  extended  some  earlier  calculations  of 
Sudan7  in  which  the  electrons  are  treated  as  a  fluid  obeying  the  following  equations: 


V2<p  =  4ir  ne,  (31) 

j  =  — -  V  x  B  =  ~VB  x  z  —  -  new,  (32) 

47r  4tt 

v-Vv  =  (</m)  (%-vxB/c)  -  Vu,  (33) 

where  U  is  the  electron  enthalpy  per  unit  mass  and  the  electron  fluid  is  assumed  to 
be  isentropic.  We  take  d/dz  =  0  and  B  =  B(x ,  y)z.  From  Eqns.  (33)  and  (32)  it 
is  easy  to  derive 

V  =  (z  x  V0)/(U  -  u),  (34) 

j  =  (V^>  x  z)ne/(fl  —  ur)  =  (c/4n)VB  x  z ,  (35) 

with  xp  =  ttp/M  —  v2/2  —  U,  fl  =  eB/mc  and  a ;  =  z  •  V  x  v  =  V  •  v  x  i  is  the 
electron  vorticity.  We  now  specialize  to  a  class  of  solutions  such  that 


n(x,y)/B(x,y)  =  no/Bo  =  constant, 


(36) 


w(x,y)/U(x,y)  =  uo/Uq  =  constant. 


(37) 
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From  (35),  (36),  and  (37)  we  conclude  that 


B  =  4nm(no/Bo)  (1  -  wo/flo)'1^, 
d2b  d2b 


dx2  +  dy 2 

d2$  d2$ 

dx 2  dy 2 


=  e 


—6 


l%-‘, 

U)Q 


(38) 

(39) 

(40) 


where  e-4  =  B/Bo,^  =  t<pjmc 2,  flou»o  <  0>  and  the  coordinates  2  and  y  have 
been  normalized  by  the  length  l  =  [Bo/(4?r  n0m)i]/|won0|5.  Equation  (39)  is  the 
Liouville  equation  with  the  solution8 


e"4  =  8[(dg2/dx)2  +  (dg2/dy)2]/(l  +  g2  +  h2)2,  (41) 

/(f)  —  fix  +  iy)  =  d{z,y)  +  ih(x,y)  is  an  arbitrary  single  valued  analytic  function 
chosen  to  fit  the  boundary  conditions.  Once  b  is  solved  the  electrostatic  potential 
may  be  obtained  from  Eqn.  (40).  The  characteristic  width  of  the  current  layer  in 
the  vacuum  gap  is 

l  =  [Bo/(47rn0m)^]/(u;ono)^  =  a~^(c/ut)  =  \e/<x$ 


where  a  =  |aJo/flo|  <  1.  Equation  (39)  is  solved  by  demanding  B  =  Bo  along  a 
profile  which  ccnnects  the  electron  sheath  surface  of  the  incoming  electromagnetic 
pulse  to  the  plasma  surface,  and  B  =  Ba(x )  along  the  plasma  surface;  B(x )  —+  0  as 
x  »  l. 

m.G  Discussion  and  Conclusions 

Our  principal  conclusions  are  as  follows:  (1)  The  renormalized  estimate  of  the 
electron  collision  frequency  is  much  less  than  ue  and  the  electron  heating  rate  is 
correspondingly  reduced.  (2)  The  snowplow  action  of  the  magnetic  field  increases 
the  plasma  density  rapidly  in  a  time  ~  2/v/0en,-  or  t2  ~  2i'«/Oen,-  causing 
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rapid  ion  heating.  The  counterstreaming  ions  lead  to  ion  viscosity  and  to  further 
ion  heating  due  to  collective  effects.  (3)  The  snowplow  velocity  is  approximately 
Bo/y/8imoM.  Near  the  cathode  and  anode  the  density  no  is  expected  to  be  less  than 
in  the  bulk  plasma.  Thus  the  snowplow  action  is  more  rapid  near  the  electrodes 
and  creates  vacuum  gaps  between  the  plasma  and  the  electrodes.  Subsequently 
snowplow  action  pushes  the  plasma  away  from  the  electrodes,  as  well  as  reducing 
the  x-dimension.  (4)  Electrons  can  cross  the  vacuum  gap  in  a  strong  magnetic 
field  by  virtue  of  a  two-dimensional  equilibrium.  When  the  electrode  gaps  become 
very  large  as  a  result  of  snowplow  action  the  boundary  conditions  for  even  a  two- 
dimensional  flow  may  not  be  satisfied  causing  current  interruption.  (5)  Plasma 
surfaces  may  be  rippled  by  the  Rayleigh- Taylor  instability  with  a  typical  wavelength 
~  current  sheath  thickness. 
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I3I.I  Appendix:  Renormalization  of  Electron  and  Ion  Equations 
for  Microinstabilities 

We  examine  nonlinear  effects  on  the  effective  electron  collision  frequency  v, 
due  to  ion  acoustic  instability  computed  by  Kulsrud  et  al.5  The  stochastic  electron 
ExB  motion  decorrelates  the  electrons  from  the  waves,  thereby  reducing  the  growth 
rate.  The  ions  are  also  subject  to  random  motion  and  the  resulting  eddy  viscosity 
.  absorbs  wave  momentum  and  saturates  the  acoustic  wave  intensity  (e^/T1)2  to  much 
less  than  unity.  The  reduced  rate  of  momentum  transfer  from  electrons  to  waves 
results  in  a  lower  i/„  than  the  linear  estimate  of  Kulsrud  et.  al.  6 

The  electrons  are  governed  by  the  drift-kinetic  equation  (|u  —  fc||V||  —  kvw\  C 

^  fle) 


df  /dt  +  \E  ■  Vf  +  V||  Vj|/  -  (e/m)E||d//di>||  =  0  (Al) 

and  the  ions  by  the  fluid  equations  neglecting  pressure.  Quasineutrality  and  the 
electrostatic  approximation  apply.  The  renormalized  drift  equation  for  the  fluctua¬ 
tions  is 


-t(u>  -  /t||t>||  -  kvw  -  kxua)f  +  i(e/m)ki\[dfo/dvn)4> 

+  kx2D.(f-f0e4>/T)  +  ^-Dv-^-f  =  0  (A2) 

dv\\  dv\\ 

where  fc2  =  fc2  4-  k2  >>  fc2,  w  and  ua  are  the  drifts  due  to  the  electrostatic  field 
given  by  Eqn.  (23)  and  by  the  drag  force  on  electrons  by  the  emitted  acoustic  waves 
respectively,  Dt  ~<j  6\e  !2  T«  >  is  the  diffusivity  of  the  nonadiabatic  electron 
fluctuations  due  to  random  ExB  motion,  and  Dv  ~<|  eE^/m  |2  re  >  is  the 
diffusivity  in  V||-space  due  to  random  acceleration  by  E\\.  Here  t~x  is  the  largest 
of  the  decorrelation  rates  due  to  (i)  electron  transit  frequency  in  the  wavepacket 
r~x  ~  &kvw,  (ii)  nonlinear  decorrelation  frequencies  k2Dt  amd  Dv/v2.  For  a 
high  level  of  fluctuations,  the  perpendicular  diffusion  is  dominant.  Thus  r"1  ~ 
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k2D,  ~  (*3v’/ne){|e^/r|)  >  Akvxv.  It  turns  out  that  this  demands  fct/e/fle  » 
w/ve.  For  these  large  decorrelation  rates  Eqn.  (A2)  may  be  approximated  by 
f  fs  (l  +  foQ-  The  electrons  behave  like  a  fluid  and  integrating  over  vj| 

we  get 

h  (  .u  —  kyw\  e<j>  . 

<A3> 

The  ions  are  unmagnetized  since  u  »  H,-  and  in  the  renormalized  ion  equations,  the 
nonlinear  term  u*  V  u  is  replaced  by  the  eddy  viscosity  — jD,V2u,  where  £>,-  ~ 

The  momentum  equation  gives  u  =  c*2  from  the  continuity  equation 

we  obtain 

*4> 


h  ~  2  k2 

no  *  u>(u>  +  ik2D{ )  T 


(A4) 


where  D,-  ~  ^y(|^r |2).  Quasineutrality  then  furnishes  the  dispersion  relation 


or 


-i 


(A5) 


k2c2  ~  V'  '  *  w  /  V*  '  '  k2De 

from  which  we  obtain  the  frequency  wr  «  fee,  and  the  renormalized  growth  rate  7 


7/wr  « 


k2Di  kyw  —  ujr 


2a >r 


+ 


2k2D„ 


(A6). 


For  comparison,  the  linear  growth  rate  7  derived  by  Kulsrud  et.  al.6  (but  valid 
only  for  we  C  fle)  is  7/fcc,  w  1.5  •  Not  only  does  the  turbulence  reduce 

the  electron  energy  input,  but  it  also  damps  the  waves  by  allowing  transfer  of 
momentum  from  waves  to  ions.  As  a  first  estimate  of  saturation  level,  set  7  =  0 
to  obtain  (j^|)  ~  (w/ve)$  (tle/kve)$  <C  1.  The  momentum  transfer  rate  from 
electrons  to  waves  is  proportional  to  the  second  term  on  the  RHS  of  (A6).  Hence, 

v-ikuUrliZ'L.)  =  3  (WVe)*(ne/&Ve)$  «  1  sinCe  M  noted  aboVe  that 

1  >  kve/Ue  >>  w/ve. 
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IV.  RESULTS  OF  EXPLORATORY  EXPERIMENTS 

IV.A  Introduction 

In  order  to  determine  the  feasibility  of  using  LION  for  gas  puff  implosions  and 
with  fine  wire  loads,  and  to  test  the  relative  soft  x-ray  yield  of  different  z-pinch 
load  configurations  with  a  constant  input  power  pulse,  two  exploratory  series  of 
experiments  have  been  performed.  Approximately  100  test  pulses  were  carried  out 
all  told.  Most  of  them  were  gas  puff  implosions,  a  few  of  which  were  carried  out 
using  a  stabilizing  magnetic  field.  About  30  shots  used  carbon  and/or  aluminum 
fiber  loads  to  form  the  z-pinch,  including  single  fiber,  x-pinch  and  “double-cross" 
configurations.  There  were  also  some  short  circuit  calibration  pulses  and  one  “open 
circuit"  test. 

The  “project-oriented”  goal  of  the  first  set  of  the  experiments  was  to  learn  to 
set  up  and  carry  out  well-diagnosed  z-pinch  experiments  with  a  full  power  pulse  on 
LION.  We  believe  we  achieved  that  project-oriented  goal:  we  obtained  useful  data 
on  45  of  70  pulses  with  a  full  set  of  diagnostics,  and  many  of  those  pulses  were  carried 
out  at  or  above  90%  of  full  generator  power.  About  40  pulses  utilized  acetylene  gas 
(C2H2)  in  the  puff  valve,  while  the  remaining  20  were  done  with  carbon  monoxide 
(CO).  About  10  tests  were  performed  with  very  fine  (6  fim  mean  diameter)  carbon 
filaments.  Included  in  the  diagnostic  package  in  the  first  experimental  series  was 
a  “McDonald”  soft  x-ray  spectrograph  which  we  borrowed  from  Sandia  National 
Laboratories,  Albuquerque  (SNLA). 

In  the  second  (25  pulse)  series  of  experiments,  the  McDonald  spectrograph  was 
not  available.  Therefore,  the  emphasis  was  on  obtaining  x-ray  yield  information  for 
different  z-pinch  configurations  driven  by  the  same  input  power  pulse. 

In  the  following  paragraphs,  we  briefly  describe  the  apparatus  and  diagnostics, 
present  detailed  results  from  one  gas  puff  test,  and  then  present  a  comparison  of 
two  fiber  pinch  tests,  one  a  “normal”  x-pinch  and  the  other  a  “double-cross.” 
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IV.B  Apparatus 


A  schematic  diagram  of  the  experiment  for  gas  puff  implosions  is  shown  in 
Figure  IV. 1.  A  puff  valve,  the  design  of  which  was  obtained  from  PI,  has  an  annular 
nozzle,  thereby  producing  a  hollow  gas  shell  with  1.5  cm  mean  radius,  and  0.5  cm 
thickness.  The  nozzle  mach  number  is  5  according  to  an  isentropic  flow  model.  The 
mass  load  in  the  gas  shell,  when  optimized  to  the  LION  power  pulse  experimentally, 
was  about  4/ig/cm  or  1017  atoms  per  cm.  This  was  estimated  from  the  puff  pressure 
measured  using  a  fast  ionization  gauge.  The  gas  from  the  supersonic  nozzle  was 
preionized  by  a  discharge  and  imploded  by  the  30  ns  rise  time,  550  kA  peak  current 
pulse  from  LION.  The  gap  between  the  anode  and  cathode  was  1.5  cm  and  the 
return  current  rods  were  4.0  cm  away  from  the  pinch  axis.  On  some  pulses,  a 
plasma  opening  switch  was  used  to  improve  the  current  rise  time. 

Experimental  diagnostics  used  in  the  first  series  of  experiments,  all  of  which 
viewed  the  plasma  from  the  side,  included:  a  streak  camera  for  time  dependent 
measurements  of  visible  light  from  the  z-pinch  or  from  x-rays  impinging  upon  a 
scintillator;  filtered  x-ray  diodes  and  a  pin  diode  to  record  x-ray  and  XUV  emission 
in  selected  spectral  regions;  a  0.5-m  spectrometer  to  measure  visible  and  UV  spectra; 
and  a  grazing  incidence  grating  spectrograph  to  collect  soft  x-ray  spectra.  (We  are 
very  much  indebted  to  Keith  Matzen’s  group  at  SNLA  for  lending  us  the  x-ray 
spectrograph  used  for  these  experiments.)  In  the  second  series  of  experiments,  the 
soft  x-ray  spectrograph  was  not  available.  However,  a  visible  light  framing  camera 
was  used. 

In  the  gas  puff  pulses,  the  amount  of  implosion  mass  was  “optimized”  by  ad¬ 
justing  the  puff  valve  plenum  pressure  such  that  the  pin-diode  showed  the  strongest 
x-ray  emission  at  or  near  the  pinch  time.  An  early  peak  on  some  pin-diode  signals 
resulted  from  hard  x-rays  early  in  the  pulse,  when  the  voltage  was  as  high  as  800 
kV. 


26 


Figure  IV. Is  Schematic  diagram  of  the  z-pinch  apparatus  when  set  up  for  gas  puff  pulses. 

There  are  eight  symmetrically  placed  return  current  rods  and  eight  side- 
on  diagnostic  ports  which  view  between  the  rods.  Preionization  of  the  gas 
puff  in  the  main  z-pinch  gap  is  accomplished  by  a  capacitor  discharged 
through  the  gas  between  the  downstream  diode  grid  and  the  preionizer 
grid. 


27 


In  order  to  exemplify  the  data  we  obtained,  we  present  details  of  a  gas  puff 
implosion  using  CO  gas,  after  which  we  will  summarize  the  key  features  of  the 
ensemble  of  data  obtained  in  the  first  set  of  experiments. 

Figure  IV.2a  shows  the  current  pulse  and  Fig.  IV.2b  the  voltage  pulse  for  the 
CO  shot.  The  pinch  occurred  60  ns  after  the  start  of  the  current  as  indicated  by  the 
changing  sign  of  the  voltage.  The  filtered  (6//m  mylar)  pin-diode  signal  is  shown  in 
Fig.  IV.3.  Note  the  half  maximum  width  of  15ns.  Figures  IV.4a  and  IV.4b  are  the 
soft  x-ray  emission  measured  by  x-ray  (XRDs)  diodes  with  0.5  and  1.0// m  aluminum 
foil  filters,  respectively.  The  transmission  window  of  these  filters  is  shown  in  Fig. 
IV.5. 

The  estimated  electron  density  and  temperature  for  the  specific  CO  implosion 
being  presented  are  of  the  order  of  1019/cm3  and  150  eV,  respectively.  These  were 
obtained  using  spectroscopy  in  the  following  way.  Figures  IV.6-IV.8  show  data 
obtained  with  the  0.5m  spectrometer  at  5292 A  (a  CVI  line),  2277A  (a  CV  line) 
and  5270A  (background),  respectively.  (The  first  and  third  wavelength  data  were 
obtained  on  the  one  shot  under  discussion  here,  and  the  second  wavelength  data 
was  obtained  on  another  virtually  identical  pulse.)  The  first  emission  peaks  shown 
in  Figs.  IV .6  and  IV.7  correspond  to  the  ionization  of  carbon  ions  through  four, 
and  then  five  times  ionized  at  high  electron  temperature  before  the  pinch  is  fully 
formed;  the  second  peak  results  from  recombination  after  the  pinch.  For  this  to 
occur  requires  the  plasma  temperature  to  be  in  excess  of  60eV. 

Strong  emission  of  resonance  line  radiation  from  H-like  and  He-like  C  and  O 
ions  was  observed  using  the  soft  x-ray  spectrometer.  Figures  IV .9  and  IV.10  show 
portions  of  time  integrated  x-ray  spectra  obtained  with  a  microdensitometer  from 
the  x-ray  film.  Because  of  the  highly  uncertain  x-ray  film  (Kodak  type  101)  calibra¬ 
tion,  electron  temperature  information  could  only  be  obtained  from  pairs  of  lines 
very  close  together  in  wavelength.  The  line  ratio  CV  34.97A/CVI  33.73A  yielded 
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an  electron  temperature  of  at  least  lOOeV,  and  the  line  ratio  of  OVIII  18.97A/ O  VII 
18.63A  yielded  at  least  140eV.  These  estimated  electron  temperatures  are  lower 
limits  because  they  correspond  to  values  if  the  plasma  were  in  equilibrium.  Thus, 
we  estimate  150eV. 

Electron  temperature  or  density  can  also  be  derived  from  the  line  ratio  of 
CV  40.73A/  CV  40.27A.  However,  because  these  are  C+4  lines,  the  density  and 
temperature  are  lower  than  those  at  pinch  time,  as  strong  radiation  in  these  lines 
occurred  before  the  pinch  time. 


Figure  IV.2a:  Current  (a.)  and  “corrected”  voltage  (b.)  waveforms  for  CO  gas  puff 
implosion  pulse  number  313.  The  negative  fiducial  mark  at  the  beginning 
of  the  voltage  waveform  corresponds  to  the  time  of  the  negative  portion 
of  the  fiducial  mark  on  the  current  waveform. 
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Figure  IV .3:  Pin  diode  trace  for  pulse  number  313.  A  6/zm  mylar  “filter”  was  used. 
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Figure  IV. 4: 


XRD  signals  for  shot  number  313.  The  XRD  result  in  (a.)  was  obtained 
with  a  0.5Mm  aluminum  filter  while  (b.)  was  obtained  with  1.0/im  of 
aluminum. 
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Figure  IV.5:  XRD  response  with  no  filter  and  with  0.5,  1.0,  and  1.5  /xm  of  aluminum 
foil  filters.  Notice  that  the  horizontal  scale  means  1.0  corresponds  to  100 
eV,  etc.  In  the  vertical  scale,  A/MW  —  Amperes/Megawatt,  is  equivalent 
to  10~6  electrons/incident  photon. 


Figure  IV.6:  Oscillograph  of  the  5292A  CVI  (C+s)  line  for  pulse  number  313. 


Figure  IV. 7:  Oscillograph  of  the  2277A  CV  (C+4)  line  from  a  shot  virtually  identical 

to  pulse  number  313. 


Figure  IV .8:  Continuum  light  at  5270A  on  pulse  number  313.  This  data  was  taken 

together  with  the  5292A  line  radiation  shown  in  Figure  7  by  the  use  of 
separate  fiber  optics  located  at  the  two  spectral  positions  in  the  focal 
plane  of  the  0.5m  spectrometer. 
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CO  Gas  Puff  Plasma 


t: 


S 

n 

o 


Wavelength  (A)  — 


Figure  TV .9:  Short  wavelength  data  from  the  x-ray  spectrograph  from  a  CO  puff  pulse. 
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Figure  IV .10:  Long  wavelength  data  from  the  x-ray  spectrometer  for  a  C2H2  puff  pulse. 

Only  the  182A  CVI  carbon  line  is  indicated  although  many  other  peaks  in 
the  spectrum  are  identified  as  carbon  lines.  In  CO  puff  pulses,  this  region 
of  the  spectrum  contains  numerous  oxygen  lines  also. 
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Knowing  the  implosion  velocity  and  the  time  delay  between  CV  line  emission 
and  the  pinch  reaching  minimum  radius,  the  estimated  pinch  density  is  about  20 
times  the  density  from  the  line  ratio  technique,  namely  about  10ie/cm3.  In  fact, 
using  the  ratio  of  intensities  of  the  equivalent  OVII  (0+6)  lines,  at  21. 6A  and 
21.8A,  which  are  radiated  much  closer  to  pinch  time,  we  do  obtain  a  density  over 
10 19 /cm3.  Figure  IV.ll  shows  a  streak  photograph  from  a  different  gas  puff  pulse. 
The  pinch  duration  inferred  from  this  photograph  is  about  10  ns.  A  2  x  107  cm/s 
implosion  velocity  is  also  implied  by  this  streak.  A  pinch  radius  of  1  mm  or  less 
is  inferred  from  streaks  using  Pilot  B  scintillator  to  convert  x-rays  to  visible  light 
which  is  then  viewed  by  the  streak  camera. 

We  now  turn  briefly  to  pulses  with  an  applied  magnetic  field,  and  then  x-pinch 
loads.  When  a  stabilizing  magnetic  field  of  3  kG  was  applied  using  the  coils  shown 
in  Fig.  IV.  1,  the  x-ray  intensity  on  the  pin-diode  was  reduced  by  a  factor  of  2.5, 
and  the  33.7 A  CVI  radiation  on  the  x-ray  film  was  estimated  to  be  a  factor  of  2 
lower.  With  a  magnetic  field  of  5  kG,  the  soft  x-ray  signals  were  very  much  weaker. 

In  the  first  experimental  series,  fiber  pinch  tests  utilized  only  carbon  fiber 
loads.  These  were  configured  in  the  form  of  an  X,  touching  in  the  middle.  As 
others  have  observed  before,1’2  extremely  bright  x-ray  sources  were  observed  at  the 
fiber  crossing  point.  The  point  x-ray  source  (about  1  mm  in  radius)  produced  about 
the  same  energy  in  line  radiation  at  CVI  33. 7A  as  the  gas  implosion  shots.  The 
highest  peak  emission  detected  at  the  pinch  time  by  the  x-ray  diodes  and  the  pin 
diode,  and  the  shortest  pinch  time  (less  than  5  ns),  were  obtained  with  this  load 
configuration.  This  suggests  that  the  x-pinch  plasma  is  more  dense  and  hotter  than 
the  gas  puff  z-pinch.  These  results  led  to  the  second  series  of  experiments,  results 
from  which  are  described  next. 
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Figure  IV.ll:  Visible  light  streak  photograph  of  a  gas  puff  implosion  using  C2H2  gas. 

Time  progresses  downward.  The  duration  of  the  “waist”  in  the  streak  is 
10ns.  The  width  of  the  waist  is  8mm.  (The  minimum  radius  observed 
when  x-rays  are  converted  to  light  with  a  fast  scintillator  and  then  viewed 
with  the  streak  camera  is  about  1mm.) 


In  the  second  series,  normal  x-pinches  and  “double-cross”  fiber  pinches  were 
tested.  The  gap  between  the  electrodes  in  which  the  fibers  were  mounted  was 
about  4cm.  To  increase  the  x-ray  radiation  yield,  combinations  of  C  (10/zm)  and 
Al  (20/zm)  fibers  were  used.  The  estimated  soft  x-ray  radiation  yield  was  several 
hundred  joules,  but  we  do  not  have  detailed  spectral  information. 

Experimental  diagnostics  included:  filtered  x-ray  diodes  (XRDs)  to  record  x- 
ray  emission  in  selected  spectral  regions;  a  visible  light  framing  camera;  and  a 
time  integrated  x-ray  pinhole  camera.  Figures  IV.  12  and  IV.  13  show  data  from 
single  and  double-cross  (C  and  Al)  fiber  shots  at  the  same  LION  Marx  generator 
charging  voltage  (80  kV).  In  these  figures,  (a)  shows  the  current  waveform  which 
peaks  at  about  500  kA  for  both  cases;  (b)  shows  the  inductively  corrected  load 
voltage,  which  is  about  900  kV  peak,  again  for  both;  (c)  shows  the  x-ray  radiation 
(>  IkeV)  measured  by  a  pin-diode  with  a  12/zm  Kimfol  filter.  The  x-ray  emission 
is  increased  by  a  factor  of  about  2  for  the  double-cross  pulse;  (d)  shows  the  soft 
x-ray  emission  (>  500eV)  detected  by  an  XRD  with  a  1/zm  Al  foil  filter,  which 
increased  about  30%;  (e)  shows  the  x-ray  radiation  detected  by  a  scintillator  with 
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7/j.m  Al  foil  filter  and  (f)  is  the  monitor  of  the  framing  camera  with  a  framing  time 
of  10ns.  The  initial  peak  on  the  pin  diode  traces  in  both  tests  (and  perhaps  on 
the  XRD  traces  as  well)  is  believed  to  be  hard  (bremsstrahlung)  radiation  from 
accelerated  electrons  when  the  gap  voltage  is  high.  Figures  IV.14  and  IV.15  show 
the  visible  framing  and  time  integrated  x-ray  images  of  the  single  and  double-cross 
fiber  shots,  respectively,  shown  in  Figs.  IV. 12  and  IV.13.  Hot  spots  at  cross  points 
due  to  harder  x-ray  emission  (~lkeV)  are  clearly  visible  when  a  7/mn  Al  filter  is 
used,  while  the  softer  x-ray  radiation  source  (>  500eV)  is  imaged  when  the  1/im 
Al  filter  is  used.  The  harder  x-ray  radiation  is  from  a  localized  region,  as  shown 
in  Figs.  IV.14  and  IV.15.  The  minimum  hot  spot  size  indicated  by  the  film  is  less 
than  300/xm,  which  is  limited  by  the  200/xm  pinhole  size.  These  data  suggest  that 
most  of  the  harder  radiation  is  from  the  cross  point (s). 
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Waveforms  from  an  x-pinch  puls,  showing  (a)  diode  current,  (b)  load  volt¬ 
age,  (c)  pm  d, ode  with  a  12nm  Kimfol  filter,  (d)  x-ray  diode  with  a  l«n 
aluminum  fi  ter,  (e)  a  Pilot  B  scintillator  with  a  7„m  A1  filter  viewed  by 
a  photomultiplier,  and  (f)  visible  light  framing  camera  monitor. 
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Shot  No.  spn459 
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Figure  IV.13:  Waveforms  for  a  “double-cross”  pulse  showing  (a)  diode  current,  (b)  load 
voltage,  (c)  pin  diode  with  a  12/un  Kimfol  filter,  (d)  x-ray  diode  with  a 
l^m  aluminum  filter,  (e)  a  Pilot  B  scintillator  with  a  7/un  Al  filter  viewed 
by  a  photomultiplier,  and  (f)  visible  light  framing  camera  monitor. 
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Figure  IV.14:  Single  cross  data  from  the  pulse  presented  in  Fig.  IV.  12,  showing  a  10ns 
visible  light  framing  photograph  (upper  image),  a  time  integrated  x-ray  pin 
hole  photograph  obtained  with  a  7/rm  Al  foil  filter  and  a  time  integrated 
x-ray  pin  hole  photograph  obtained  with  a  l,um  Al  foil  filter. 


39 


r 


Double  Cross  Al  C  Fibers 
Visible  Framing  Camera 


X-RAY 

0.5mm  pin  hole; 
7pm  Al  foil; 

Two  hot  spots  are 
10mm  apart 


X-RAY 

0.2mm  pin  hole; 
lpm  Al  foil; 


Hot  spot  sue  <  300pm 


/  > 


Figure  IV. 15:  Double  cross  data  from  the  pulse  presented  in  Fig.  IV.  13  showing  a  10ns 

visible  light  framing  photograph  (upper  image),  a  time  integrated  x-ray  pin 
hole  photograph  obtained  with  a  7/xm  Al  foil  filter  and  a  time  integrated 
x-ray  pin  hole  photograph  obtained  with  a  lfj.m  Al  foil  filter.  For  scale, 
the  two  hot  spots  in  the  middle  and  lower  photographs  are  10mm  apart. 
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